A test for cosmic distance duality 
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Testing the cosmic distance duality relation (CDDR) constitutes an important task for cosmology 
and fundamental physics since any violation of it would be a clear evidence of new physics. In 
this paper, we propose a new test for the CDDR using only current measurements of the gas mass 
fraction of galaxy clusters from Sunyaev-Zeldovich (/sze) and X-ray surface brightness (fx-ray) 
observations. We show that the relation between fx-ray and fszE observations is given by fszE = 
nf x-ray, where n quantifies deviations from the CDDR. Since this latter expression is valid for the 
same object in a given galaxy cluster sample, the method proposed removes possible contaminations 
from different systematics error sources and redshift differences involved in luminosity and angular 
diameter distance measurements. We apply this cosmological model-independent methodology to 
the most recent fx-ray and fszE data and show that no significant violation of the CDDR is found. 

PACS numbers: 98.80.-k, 98.80.Es, 98.65.Cw 



I. INTRODUCTION 

The so-called cosmic distance duality relation 
(CDDR) pj, which is closely connected with the Ether- 
ington reciprocity theorem 0, Q > plays an important role 
in observational cosmology ranging from gravitational 
lensing studies, galaxy and galaxy clusters observations, 
to analyses of the cosmic microwave blackbody radiation 
(CMB)[4]. It relates the luminosity distance Dr, with 
the angular diameter distance Da through the identity 



£< 1 + =>- 3 



(1) 



The above result can be easily demonstrated for usual 
Friedmann-Lemaitre-Robertson- Walker cosmologies. As 
is well known, it only requires source and observer to 
be connected by null geodesies in a Riemannian space- 
time and cosmological conservation of the number of pho- 
tons p], [|| . Examples of non-standard scenarios that vio- 
late the equality ([T]) are models in which photons do not 
travel on unique null geodesies, models with variations of 
fundamental constants or with photon non-conservation 
due to coupling to particles beyond the standard model 
of particle physics, absorption by dust, etc. (see, e.g., 
[9l-frij for a discussion). 

Attempts to actually test the CDDR validity from 
astronomical observations have started only recently 
[EBUl- For instance, in Ref. @ it was argued that the 
Sunyaev-Zel'dovich effect (SZE) plus X-ray technique for 
measuring galaxy cluster distances is strongly dependent 
on the CDDR validity. As is well known, one can consider 
different electronic density dependencies combined with 
some assumptions about the galaxy cluster morphology 
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in order to evaluate its angular diameter distance with 
basis on Eq. (JTJ, i.e. 
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where D A ata (z) is the angular diameter distance to a 
given galaxy cluster obtained directly from a jointly anal- 
ysis of its SZE and X-ray surface brightness observations, 
Sxo is the central X-ray surface brightness, T e o is the 
central temperature of the intra-cluster medium, A e no 
is the central X-ray cooling function of the intra-cluster 
medium, ATo is the central decrement temperature, and 
9 C refers to a characteristic scale of the cluster along the 
line of sight (l.o.s.), whose exact meaning depends on 
the assumptions adopted to describe the galaxy cluster 
morphology. 

However, if one assumes a more general expression for 
CDDR, such that, 
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it is possible to show that (for details see (El H2j|) 

D^{z)=D A {z)r 1 {zf. 
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Therefore. 



D A ata (z) reduces to the real angular diame- 



ter distance only when the DD relation is strictly valid 
(r) = 1). In order to quantify the r\ parameter, the au- 
thors in Ref. 01 fixed Da{z) by using the cosmic concor- 
dance model [l6| while for D^ ata (z) they considered the 
18 galaxy clusters from the Ref. pj| for which a spheri- 
cally symmetric cluster geometry has been assumed. By 
assuming r\ constant, their statistical analysis provided 
r\ = 0.9llg'o4 (lc), and is therefore only marginally con- 
sistent with the standard result. 

On the other hand, the CDDR should be tested only 
from astronomical observations, i.e., finding cosmological 
sources whose intrinsic luminosities and intrinsic sizes are 
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known. Thus, after measuring the source redshift, one 
can determine both Dl and Da to test directly the rela- 
tion. In principle, this ideal method should not use any 
relationship coming from a specific cosmological model, 
i.e., they must be determined by means of intrinsic as- 
trophysical quantities only. 

In recent papers, the validity of the CDDR has been 
discussed using Da measurements from galaxy clusters 
(GC) and luminosity distances from type la supernovae 
(SNe la) @. In such analyses, subsamples of GC and 
SNe la are built so that differences in redshift between 
objects in each sample are small (Az ~ 10 -3 ), thereby 
allowing a validity test of the CDDR. There are, however, 
at least three aspects that should be considered when per- 
forming this kind of analysis. First, that estimates of 77 
from this method are potentially contaminated by differ- 
ent systematics error sources in GC and SNe la observa- 
tions. Second, that some SNe la light curve fitters use a 
specific cosmological scenario in their calibration process, 
which makes the GC/SNe la test not completely model- 
independent (see Ref. [l2j for a discussion on influence 
of the SNe la light curve fitters on the CDDR test). Fi- 
nally, that the use of different objects to derive Dl and 
Da implies necessarily in a choice for Az, which affects 
the resulting estimates of the 77 parameter (see [l3| for a 
recent discussion. See also [l4[ for an interesting model- 
independent CDDR test based on future observations of 
a redshifted 21cm signal from disk galaxies). 

In order to circumvent these observational problems, 
we propose a consistent model-independent test for Eq. 
(H|) that uses only current observations of the gas mass 
fraction of GC's. To perform our analysis, we use a sam- 
ple of fx~ray and fszE measurements of 38 GC's as 
discussed in Ref. [13. We show that when the CDDR 
is taken into account current fx -ray and fszE measure- 
ments provide a direct test for the CDDR. 



II. GAS MASS FRACTION 
A. X-ray Observations 

The gas mass fraction is defined as / — Mgas / Mxot 
[l9l |. where Mxot is the total mass obtained via hydro- 
static equilibrium assumption and M gas (gas mass) is 
obtained by integrating the gas density model. The gas 
density model frequently used is the 3-dimcnsional elec- 
tron number density of the spherical /3 model given by 

M 1 



where n e is the electron number density, r is the radius 
from the center of the cluster, r c is the core radius of the 
intracluster medium (ICM) , and /3 is a power law index. 
Under this assumption, M to t and M gas are given by (2i| 



Af tota i(r) 



G[ira p r\ 



(6) 



where T c is the temperature of the intra cluster medium, 
fx and m p are, respectively, the total mean molecular 
weight and the proton mass, fee the Boltzmann constant, 
and 



M gas {r) = A 
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where A = fi e n eo m p D A 3 , 9 = t/Da and /x e is the 
mean molecular weight of the electrons. 

The X-ray surface brightness over some frequency band 
is written as 



S x = 
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where the integral is along the line of sight and A eH is the 
X-ray cooling function, proportional to Tj^ 2 [22J. Note 
that the above equation depends explicitly on the ratio 
between the angular diameter and luminosity distances. 
Thus, if the CDDR is taken as valid [77 = 1 in Eq. ([T])] and 
a fiducial angular diameter distance, D* A , is assumed, the 
central electron density n e o can be analytically obtained, 
i.e., FH 
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providing the well-known relation in current X-ray gas 
mass fraction measurements f x-ray oc D* A 3 ^ 2 [lTj ] (in the 
above expression S x q is the central surface brightness and 
r(x) is the Gamma function). 

However, if the validity of the CDDR is not previously 
assumed (77 ^ 1), Eq. (0) is now rewritten as 
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which clearly shows that gas mass fraction measurements 
extracted from X-ray data are affected by a possible vi- 
olation of the CDDR scaling as 



f X-ray « VD A 



3/2 
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B. Sunyaev-Zel'dovich Observations 



1 It is worth mentioning that we have used f gaa data for the non- 
isothermal double /3-model in our statistical analyses. The spher- 
ical /9 model is used in this section only for simplicity but without 
loss of generality for the test proposed. 



The measured temperature decrement ATsze of the 
CMB due to the Sunyaev-Zel'dovich effect [23| is given 
by[H 



ATp 

7bMB 



/(i/,T e )^§ / n e T c dl, (12) 
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FIG. 1: a) Gas mass fraction measurements as a function of redshift for 38 galaxy clusters [TfJ. Blue circles (red squares) stand 
for measurements obtained via current X-ray (SZE) observations. Open circles and squares correspond to galaxy clusters that 
present large reduced \ 2 (Xd.o.f. > 2) when described by the hydrostatic equilibrium model, b) The likelihood distribution 
functions for the CDDR parameter 770. Blue solid lines correspond to the results for PI whereas red dashed lines for P2. In 
this analysis, the entire sample of Ref. [l5l ] was used (38 clusters), c) The same as in Panel lb for the reduced galaxy cluster 
sample (29 clusters). When compared with the previous analysis, the compatibility of the data with the CDDR is clearly more 
evident. 



where Tcmb = 2.728 K is the present-day temperature 
of the CMB, <tt the Thompson cross section, m c the 
electron mass and f{y,T e ) accounts for frequency shift 
and relativistic corrections [24{ . It worth mentioning that 
the gas temperature, T e , is insensitive to the validity of 
the CDDR since it is obtained through the shape of the 
X-ray spectra (bremsstrahlung spectra) and not through 
X-ray luminosity. 

Using SZE observations, the central electron density 
can now be expressed as 



l e0 



AT Q m e c 2 T(p) 



(13) 

which is insensitive to the validity of the CDDR. There- 
fore, current gas mass fraction measurements via SZE 
depend only on the angular diameter distance as 



fszE oc Da- 



(14) 



C. fsZE/fli 



relation 



Current fx -ray measurements have been obtained by 
assuming the validity of the CDDR [35], [36[. If, how- 
ever, this is not the case, the real gas mass fraction from 
X-ray observations (Eq.(ll)) should be related with the 
current observations by fx-ray ~ Vfx-ray So, if all 
the physics behind the X-ray and SZE observations are 
properly taken into account, one would expect f gas mea- 
surements from both tecniques to agree with each other 



since they are measuring the very same physical quan- 
tity. In this way, it is clear that the general expression 
relating current X-ray and SZE observations is given by: 



fsZE = fx-ray = Vfx- 



(15) 



which provides a direct test for the CDDR. Since the 
above expression holds for a given object, a possible in- 
fluence on the 77 estimates due to redshift differences of 
distinct objects (e.g., in tests involving SNe la and GC) 
is fully removed. 

Before discussing our estimates on the rj parameter 
from the above relation, it is worth mentioning that in 
both methods to test the CDDR (the one of Eq. (fl"5|) 
and those using Da measurements of galaxy clusters from 
their SZE and X-ray observations (jllll E2J) other physi- 
cal effects could be present in the SZE observations, lead- 
ing to results with 77 ^ 1 that not necessarily would be 
related to a CDDR violation in X-ray frequency. For 
instance, as mentioned earlier [see Eq. (13)], SZE ob- 
servations are independent of the CDDR validity. How- 
ever, these observations are redshift-independent only if 
there is no process of energy injection into the CMB. 
Otherwise, the standard linear relation of the CMB tem- 
perature evolution, Tqmb{z) = Tcmb(z = 0)(1 + z), is 
not valid 25[ and changes are needed in Eqs. (12) and 
(13). In this line, there are recent analyses testing the 
standard evolution of the CMB temperature through dif- 
ferent techniques by using a more general expression such 
asTcMB{z) = TcMB{z = 0){l + z) 1+ P [US]. As a ba- 
sic result, all analyses confirm the standard relation, i.e., 
/3 « (see Ref. [30] for recent constraints on /3 value by 
using SZE observations and carbon monoxide excitation 
at high-z). 



4 



In what follows, we explore Eq. (|15l) and discuss con- 
straints on the CDDR from current X-ray and SZE gas 
mass fraction measurements. 



III. GALAXY CLUSTER DATA 

To discuss the potential of Eq. (|T5|) in probing the 
CDDR, we use the most recent X-ray and SZE gas mass 
fraction measurements to date, as given in Ref. [15J . The 
sample consists of 38 massive galaxy clusters spanning 
redshifts from 0.14 up to 0.89 whose X-ray data were 
obtained from the Chandra X-ray Observatory and SZE 
data from the BIMA/OVRO SZE imaging project, which 
uses the Berkeley-Illinois-Maryland Association (BIMA) 
and Owens Valley Radio Observatory (OVRO) interfer- 
ometers to image the SZE. In order to perform a realistic 
model for the cluster gas distribution and take into ac- 
count a possible presence of cooling flow, the gas density 
was modeled with the non-isothermal double /3-modcl 
that generalizes the single /3-model profile given by Eq. 

Another important aspect worth mentioning is that 
the shape parameters of the gas density model (9 C and /3) 
were obtained from a joint analysis of the X-ray and SZE 
data [3l|, which makes the SZE gas mass fraction not in- 
dependent. However, as has been shown from current 
simulations [32|, the values of 9 C and f3 computed sepa- 
rately by SZE and X-ray observations agrees at la level 
within a radius r25oo (at which the mean enclosed mass 
density is equal to 2500 cosmological critical density), the 
same used in the La Roque et al. observations. There- 
fore, we believe that, unless a violation of the CDDR may 
affect even the shape parameters obtained from the X- 
ray data (and not only S x o), one can use the La Roque et 
al. sample to perform the method discussed earlier. In 
Figure la we show the sample of 38 measurements of the 
gas mass fraction of galaxy clusters obtained via X-ray 
surface brightness ( assuming 77 = 1) and the SZE (see 
also Table 5 of Ref. QJ]). 

IV. ANALYSIS AND RESULTS 

We modify the CDDR to test for any violation by using 
simple parameterizations for 77. In order to take into ac- 
count a possible influence of different 77 parameterizations 
on the results, we use in our analyses two functions fl8j : 

( 1 + Voz (PI) 

T]( Z ) = < 

[l + r] z/(l + z) (P2) 

PI is a continuous and smooth one-parameter linear ex- 
pansion, whereas P2 includes a possible epoch- dependent 
correction, which avoids the divergence at very high z. 
Note also that, differently from constant 77 parameteri- 
zations (e.g., D L /D A = 770(1 + zf 0), both PI and P2 



recover the CDDR in the limit z — > 0, as obtained from 
cosmographic derivations. 

We evaluate our statistical analysis by defining the like- 
lihood distribution function C oc e~ x / 2 , where 

x 2 = j^ m-mMz)\\ (16) 

%— 1 ^i.obs 

Vi,obs(z) = ! sze I ! x-ray and a\ ohs is the uncertainty as- 
sociated to this quantity. Also, the i index indicates the 
sum over each cluster. In Fig. lb we plot our first con- 
straints on the CDDR from Eq. (|15p . By considering 
all 38 galaxy clusters of the sample discussed earlier, we 
show the likelihood distribution as a function of the pa- 
rameter 770 for PI (blue solid line) and P2 (red dashed 
line). For these two cases, we find 770 = —0.15 ± 0.07 
(X 2 d . . f . = 1.02) and 770 = -0.22 ± 0.10 ( X L./. = 1-04) 
at la level, respectively. Note that the CDDR (770 = 0) 
is slightly compatible with theses data, being ~ 2a off 
from the best-fit values in both cases. For the sake of 
completeness and also to verify the effect of the gas mod- 
elling on the results, we also performed the same analysis 
using the isothermal j3 model. In this case, we obtained 
?7o = -0.08 ± 0.16 and 770 = -0.13 ± 0.21 at 2a level for 
the linear and non-linear parametrizations, respectively. 

An important aspect concerning the galaxy cluster 
sample shown in Fig. la is that some objects present 
questionable reduced x 2 (2.43 < Xd.o.f. — 41.62) when 
described by the hydrostatic equilibrium model (see Ta- 
ble 6 in Ref. [34]). They are: Abell 665, ZW 3146, 
RX J1347.5-1145, MS 1358.4 + 6245, Abell 1835, MACS 
J1423+2404, Abell 1914, Abell 2163, Abell 2204. By ex- 
cluding these objects from our sample (we end up with 
a subsample of 29 galaxy clusters), we perform a new 
analysis whose results are displayed in Fig. lc. We note 
that, when compared with the previous analysis (Fig. 
lb), the compatibility of the data with the CDDR is 
clearly more evident now, with 770 = —0.06 ± 0.07 (PI) 
and 770 = -0.07 ± 0.12 (P2) at la level. It is worth 
observing that for all analyses performed in this paper 
a negative value for the CDDR parameter was prefered 
by the data (although the data are fully consistent with 
770 = 0). A possible explanation for that has been dis- 
cussed in Ref. [l(| in terms of cosmic opacity or the exis- 
tence of axion-like and mini-charged particles (see, e.g., 
11] for a recent review on these weakly-interacting-sub- 
eV particles). In Table I we summarize the main results 
of our analyses. 



V. CONCLUSIONS 

In this Letter, we have discussed how current measure- 
ments of gas mass fraction of galaxy clusters from X- 
ray and SZE effect observations can be used to test the 
CDDR. We have shown that if this relation is consis- 
tently taken into account current fx -ray and fszE mea- 
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38 galaxy clusters data 


x'/d.o.f 


PI 7/0 = -0.15 ±0.14 
P2 7/0 = -0.22 ± 0.21 


1.02 
1.04 


29 galaxy clusters data 


X A /d.o.f 


PI rio = -0.06 ± 0.16 
P2 7/o = -0.07 ± 0.24 


0.93 
0.92 



TABLE I: Constraints on the CDDR parameter 7/0. The error 
bars correspond to 2a (Ax 2 = 4) . 



surements are related by fszE = vf x-ray [Eq. (fT5|) ]. 
allowing a direct test for the CDDR. 

To perform our analyses we have considered two dis- 
tinct forms for Tj(z), i.e., 7/ = I + t/qz and 7/ = l±r/oz/(l + 
z), which recover the equality between Dl and Da at 
very low redshifts. By considering 38 gas mass fraction 
measurements of galaxy clusters obtained from their SZE 
and X-ray emissions, we have found no significant influ- 
ence of the above 77(2) parameterizations on the results 
and that a no violation of the CDDR is compatible at 
2c7 level. However, if our sample is corrected to account 
for possible statistical error sources due to galaxy clus- 
ter modeling, a value compatible with the validity of the 



CDDR is found (see Table I). 

Finally, it is worth mentioning that, irrespective of the 
fiducial model adopted in observations, f gas measure- 
ments from both tecniques (X — ray and SZE) for a 
given cluster must agree and if they do not the detected 
difference must be associated with the duality parameter 
r/. Therefore, the method here proposed is cosmological 
model-independent and also presents a clear advantage 
over tests involving different kinds of observations (e.g., 
GC and SNe la) since the above relation is analyzed for 
the same object in a given sample. We believe that when 
applied to upcoming observational data the method dis- 
cussed here may be useful to probe a possible violation 
of the CDDR. 
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